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Abstract 

Although QCD can be treated perturbatively in the high energy Umit, lower 
energies require uses of nonperturbative methods such as ADS/CFT and/or 
Abelian reduction. These methods are not equivalent. While the first is re- 
stricted to supersymmetric Yang-Mills model with number of colors going to 
infinity, the second is not restricted by requirements of supersymmetry and 
is designed to work in physically realistic limit of finite number of colors. In 
this paper we provide arguments in favor of the Abelian reduction methods. 
This is achieved by further developing results of our recent works reanalyzing 
Veneziano and Veneziano-like amplitudes and the models associated with these 
amplitudes. It is shown, that the obtained new partition function for these 
amplitudes can be mapped exactly into that for the Polychronakos-Frahm (P- 
F) spin chain model recoverable from the Richardon-Gaudin (R-G) XXX spin 
chain model originally designed for treatments of the BCS-type superconductiv- 
ity. Because of this, it is demonstrated that the obtained mapping is compati- 
ble with the method of Abelian reduction. The R-G model is recovered from 
the asymptotic (WKB-type) solutions of the rational Knizhnik-Zamolodchikov 
(K-Z) equation. Linear independence of these solutions is controlled by de- 
terminants whose explicit form (up to a constant) coincides with Veneziano 
(or Veneziano-like) amplitudes. In the simplest case, the determinantal condi- 
tions coincide with those discovered by Kummer in 19-th century. Kummer's 
results admit physical interpretation by relating determinantal formula(s) to 
Veneziano-like amplitudes. Furthermore, these amplitudes can be interpreted 
as Poisson-Dirichlet distributions playing central role in the stochastic theory 
of random coagulation-fragmentation processes. Such an interpretation is com- 
plementary to that known for the Lund model widely used for description of 
coagulation- fragmentation processes in QCD. 
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1 Introduction 



Since times when quantum mechanics (QM) was born (in 1925-1926) two seem- 
ingly opposite approaches for description of atomic and subatomic pliysics were 
proposed respectively by Heisenberg and Schrodinger. Heisenberg's approach is 
aimed at providing an affirmative answer to the following question: Is combina- 
torics of spectra (of obsevables) provides sufficient information about micro- 
scopic system so that dynamics of such a system can be described in terms of 
known macroscopic concepts? Schrodinger's approach is exactly opposite and 
is aimed at providing an affirmative answer to the question: Using some plau- 
sible mathematical arguments is it possible to find equations which under some 
prescribed restrictions will reproduce the spectra of observables? Although it 
is widely believed that both approaches are equivalent, already Dirac in his 
lectures on quantum field theory, Ref. [1], noticed (without much elaboration) 
that Schrodinger's description of QM contains a lot of " dead wood" which can 
be safely disposed. According to Dirac "Heisenberg's picture of QM is good 
because Heisenberg's equations of motion make sense" . 

To our knowledge, Dirac's comments were completely ignored, perhaps, be- 
cause he had not provided enough evidence making Heisenberg's description of 
QM superior to that of Schrodinger's. In recent papers, Ref. [2,3], we found 
examples supporting Dirac's claims. From the point of view of combinatorics, 
there is not much difference in description of QM, quantum field theory or 
string theory as demonstrated in Ref. [4]. Therefore, in this paper we adopt 
Heisenberg's point of view at QCD and string theory using results of our recent 
works aimed at re analyzing the existing treatments connecting Veneziano (and 
Veneziano-like) amplitudes with the respective string-theoretic models. As re- 
sult, we found new tac;hyon-free models reproducing Veneziano (and Veneziano- 
like) amplitudes. In this work results of our papers, Ref.s.[5-7], to be called 
respectively as Part I, Part II and Part III are developed further. This has be- 
come possible in view of the work by Reshctikhin and Varchenko, Ref. [8], and 
by Varchenko summarized in Varchenko's MIT lecture notes, Ref. [9]. They 
enabled us to relate Veneziano (and Veneziano-like) amplitudes (e.g. those de- 
scribing TTTT scattering) to Knizhnik-Zamolodchikov (K-Z) equations and, hence, 
to WZNW models. This is achieved by employing known connections between 
WZNW models and spin chains. In the present case, between the K-Z equations 
and the XXX-type Richardson-Gaudin (R-G) spin chains. To keep things in a 
perspective, we would like to provide some rationale behind the above mentioned 
connection with spin chains developed in this paper. 

As is well known, all information in high energy physics is obtainable through 
proper interpretation of scattering data. It is believed that for sufficiently high 
energies such data are well described by the phenomenological Regge theory 
and can be conveniently summarized with help of Chew-Frauthchi (C-F) plots 
relating masses to spins (angular momenta), e.g. see book by Collins, Ref. [10]. 
Veneziano amplitudes are by design Regge-behaving. Both Regge theory and 
Veneziano amplitudes emerged before major developments in QCD took place 
in 70ies. Once these developments took place, naturally, it was of interest to 
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recover the Regge theory from QCD. Even though there are many ways of doing 
so, to our knowledge, the problem is still not solved completely. This is so for 
the following reasons. 

Although the amount of data obtained by perturbativc treatments of QCD 
is quite impressive, e.g. read Ref.[ll], these results are not helpful for estab- 
lishing the Regge-type behavior of QCD. Since such a behavior can be easily 
established with help of variety of string models, the task lies is connecting these 
models with QCD. For the sake of space, we do not discuss broader reasons (e.g. 
quantum gravity) for development of string theory. Recently, another path- 
ways towards establishing Regge-type behavior were explored. For instance, in 
Ref.[12], in the large N limit (to be defined in the main text) by ignoring masses 
of quarks the spin zero glueball mass spectrum was obtained analytically for 
2+1 dimensional pure Yang- Mills theory^ which in the high energy limit fits 
perfectly Regge theory. Although in Ref.[12] string models were not used, the 
obtained results are in excellent agreement with large N lattice calculations. 
Given this, they still suffer from several drawbacks. First, the large N limit, 
even though mathematically convenient, physically is questionable as it will be 
explained below, in the main text. Second, more physically interesting are the 
spectra of mesons and baryons. These spectra are traditionally obtained with 
help of string-theoretic models. Since calculations involving these models are 
by design made Regge-behaving, the task was (still is) to connect these string 
models with QCD. It is widely believed that such a connection is achievable 
via ADS/CFT correspondence [14] between strings living in anti-de Sitter space 
and A/" = 4, N— > oo supersymmetric Yang-Mills (Y-M) model in = 3 -|- 1 di- 
mensions. By design, such supersymmetric Y-M model does not contain quark 
masses'^. Again, such correspondence becomes physically meaningful only if 
the number of colors N could be made finite and the supersymmetry could be 
broken. Unfortunately, these requirements are in apparent contradiction with 
the way the ADS/CFT correspondence was established in the first place thus 
making such a task very difficult to accomplish. 

In 1981 't Hooft suggested in Ref.[18] to reduce the non Abelian QCD to 
Abelian Ginzburg-Landau (G-L) type theory. The rationale for such an Abelian 
reduction can be traced back to the work of Nambu, Ref. [19] . In his work Nambu 
superimposed G-L theory with the theory of Dirac monopoles to demonstrate 
quark confinement for mesons. Incidentally, Veneziano amplitudes are suited 
the most for describing meson resonances, e.g. see Ref. [10]. These are made 
of just two quarks: quark and antiquark. Thus, if the existence of Abelian 
reduction would be considered as proven, this then would be equivalent to the 
proof of quark confinement. Recent numerical studies have provided convincing 
evidence supporting the idea of quark confinement through monopole condensa- 
tion, e.g. see Refs. [20,21]. Since publication of 't Hooft's paper many theoretical 
advancements were made, most notably by Cho, Ref.s[22,23], and Kondo, Ref.s 

^Thcsc results were subsequently extended to 3-|-l dimensions in Ref. [13]. 

■^Thc achievents of this method are summarized in excellent recent reviews by Benna and 
Klebanov, Ref. [15], and Brodsky, Ref. [16]. Our work, Ref. [17], provides an introduction to 
the ideas and methods of ADS/CFT. 
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[24-26], whose work was motivated by that by Faddeev and his group. Results 
of this group are summarized in the recent review by Faddeev, Ref. [27]. From 
this reference it follows that most of efforts to date were spent on description of 
the masslcss version of QCD (just like in the case of ADS/CFT correspondence, 
but without invoking supersymmetry or requiring N— > oo). Excitation spectrum 
of solitonic knotted-like structures (admitting interpretation in terms of closed 
strings) provides the spectrum of glueball masses as demonstrated in Ref. [25]'*. 
In a recent paper, Ref. [28], Auckly, Kapitanski and Speight demonstrated how 
Skyrme model can be obtained from Faddeev model. Since Skyrme model was 
used for a long time for description of the baryon spectra, e.g. read Ref. [29], and 
since already Nambu recognized usefulness of the Abelian reduction for descrip- 
tion of meson spectra, it follows that the Abelian reduction method is capable of 
providing sufficient information about QCD in the strong coupling regime. Fur- 
thermore, in another paper, Ref. [30], Kapitanski and Auckley obtained result of 
major importance for this paper. They demonstrated that it is always possible to 
find such Chern-Simons (C-S) fimctional which upon minimization will produce 
the same results as those obtained by minimization of either Skyrme or Faddeev 
model. Asorey, Falceto and Sierra in Ref. [31] demonstrated how C-S model is 
related to microscopic BCS model of superconductivity. Since the Abelian re- 
duction of QCD produces G-L-type model (or collection of G-L models^), and 
since the underlying microscopic model is of BCS-type whose elementary exci- 
tations are described by the Richardson-Gaudin spin chain model, Ref. [31], the 
task of this work lies in demonstrating that the combinatorics of scattering pro- 
cesses associated with Veneziano (and/or Veneziano-like) amplitudes leads to 
the R,-G spin chain model of BCS superconductivity. The demonstrated in this 
work spin chain connection made in the spirit of Heisenberg's work on quantum 
mechanics, Ref. [2], favors the Abelian reduction method over ADS/CFT. Un- 
like other Schrodinger-style papers discussed above-all extracting the Abelian 
Ginzburg-Landau-type model from the non Abelian QCD- the results of this 
work use only combinatorics of scattering data as an input to arrive at the same 
conclusions. 

It should be noted that connections between either QCD and spin chains or 
between strings^ and spin chains were already discussed in literature for quite 
some time. Recent paper by Dorcy, Ref. [32], contains may references listing 
these earlier results. Subsequently, they had been replaced by those whose 
methods are based on ADS/CFT correspondence. From the point of view of 
this correspondence, connections between strings and QCD also can be made 
through spin chains as it is demonstrated in the seminal paper by Gubser, 
Klebanov and Polyakov, Ref. [33]. Their ideas were developed in great detail 
in the paper by Minahan and Zarembo, Ref. [34]. The spectrum of anomalous 
dimensions of operators in the TV = 4, N—>- oo supersymmetric Yang-Mills (Y- 

"^It should be noted that the problem of existence of a mass gap in QCD was included into 
seven the most outstanding millciiium prize mathematical problems of the 21st century by 
the Clay Mathematical Institute, e.g. see http://www.claymath.org/millenmum/ 

®E.g. read Faddeev's paper, Ref. [27]. 

6 Within the context of ADS/CFT. 
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M) model (described in terms of excitation spectrum of the spin chain model) is 
related to the string spectrum describing hadron masses. To connect these facts 
with developments in this paper we mention papers by Kruczenski, Ref.[35], and 
Cotrone et al, Ref.[36]. In both papers spin 1/2 XXX Hcisenbcrg chain was used 
for description of excitation spectrum. Furthermore, in the paper by Cotrone et 
al explicit connection with the hadron mass spectrum was made. Both papers 
invoke ADS/CFT correspondence. In this work, we reobtain these spin chain 
results using combinatorial arguments following Heisenberg's philosophy''. 

In Section 2 we reobtain Veneziano partition function (derived much more 
rigorously in Part II) . In Section 3 we demonstrate that this partition function 
coincides with that for the Polychronakos-Frahm (P-F) spin chain model. Al- 
though such a model was studied extensively in literature, in Ref.[37] we discuss 
a variety of new pathways establishing links between the P-F spin chain and 
many known string-theoretic models, including the most recent ones. For the 
sake of space, we do not reproduce these results in this work. Instead, in Section 
4 we use results of Rcshetikhin and Varchcnko, Rcf. [8], and Varchenko, Rcf. [9], 
in order to make a connection between the P-F and R-G XXX spin chains. 
We use the results of our Part II in order to demonstrate that the excitation 
spectrum of new Veneziano model obtained in Part II coincides exactly with 
that for the R-G spin chain. 

Since both ADS/CFT and this work point towards the same spin 1/2 XXX 
chain, this cannot be considered as purely coincidental. Section 5 provides 
some explanations of noticed coincidence based on independent combinatorial 
arguments having their origin in the theory of random fragmentation and coag- 
ulation processes summarized in Ref.s[38-40]. This theory was adopted for high 
energy physics by Mekjian, e.g. see Ref.[41] and references therein. A different 
approach to coagulation-fragmentation processes in QCD is developed by An- 
dersson and collaborators and is known in literature as the Lund model, Ref. [42] . 
In this work no attempts are made to compare these two approaches. Instead, 
in Section 5 we argue that in the theory of coagulation-fragmentation processes, 
Veneziano amplitudes play the central role. In this theory they are known as 
the Poisson-Diriclet (P-D) probability distributions. The discrete spectra of all 
exactly solvable quantum mechanical (QM), field and string-theoretic models 
can be rederived/reobtained in terms of the obscrvablcs for the P-D stochastic 
processes. This is so because all exactly solvable QM problems involve some 
kind of orthogonal polynomials-all derivable from the Gauss hypergeometric 
function-admiting an interpretation in terms of the P-D process. Since the K-Z 
equations are essentially the hypergeometric equations of many variables, e.g. 
see Ref. [9] , and since all nontrivial Feynman diagrams of quantum field theory 
can be looked upon as solutions of these hypergeometric equations as explained 
in Ref. [4], the processes they describe are also of the P-D type. 

Finally, Appendix contains several auxiliary results complementing those 
presented in Sections 4 and 5. 

review of Heisenberg's arguments resulted in birth of moden quantum mechanics can 
be found in Ref. [2]. 
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2 Combinatorics of Veneziano amplitudes and 
spin chains. Qualitative considerations 

In Part I, we noticed that the Veneziano condition for the 4-particle amphtude 
given by 

a{s) + a{t) + a{u) = -1, (2.1) 

where a(s), a{t), a{u) £ Z, can be rewritten in a more suggestive form. To this 
purpose, following Ref. [43] , we need to consider additional homogenous equation 
of the type 

a{s)m + a{t)n + a{u)l + fc • 1 = (2.2) 

with m,n,l,k being some integers. By adding this equation to Eq.(2.1) we 
obtain, 

a{s)'m + a{t)n + a{u)l = k (2.3a) 

or, equivalently, as 

m + n2 + ns = N, (2.3b) 

where all entries by design are nonnegative integers. For the multiparticle case 
this equation should be replaced by 

no + • • • + rife = TV (2.4) 

so that combinatorially the task lies in finding all nonnegative integer combi- 
nations of no, ...,nk producing Eq.(2.4). It should be noted that such a task 
makes sense as long as N is assigned. But the actual value of A'' is not fixed and, 
hence, can be chosen quite arbitrarily. Eq.(2.1) is a simple statement about the 
energy-momentum conservation. Although the numerical entries in this equa- 
tion can be changed as we just explained, the actual physical values can be 
subsequently reobtained by the appropriate coordinate shift. Such a procedure 
should be applied to the amplitudes of conformal field theories (CFT) with 
some caution since the periodic (or antiperiodic, etc.) boimdary conditions 
cause energy and momenta to become a quasi- energy and a quasi momenta (as 
it is known from solid state physics). 

The arbitrariness of selecting N reflects kind of a gauge freedom. As in gauge 
theories, we may try to fix the gauge by using some physical considerations. 
These include, for example, an observation made in Part I that the four particle 
amplitude is zero if any two entries into Eq.(2.1) are the same. This fact causes 
us to arrange the entries in Eq.(2.3b) in accordance with their magnitudes, e.g. 
ni>n2>n^- More generally, we can write: no > ni > • • • > > 1^. 

If the entries in this sequence of inequalities arc treated as random nonneg- 
ative numbers subject to the constraint given by Eq.(2.4), these requirements 
are necessary and sufficient for recovery of the probability density for such set 
of random numbers. This density is known in mathematics as the Dirichlet 

^The last inequality: n/, > 1, is chosen only for the sake of compajrison with the existing 
literature conventions, e.g. see Ref. [44]. 
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distribution^ [38-40,45]. Without normalization, integrals over this distribu- 
tion coincide with Veneziano amplitudes. Details are given in Section 5. Thus, 
Veneziano condition leads to Veneziano amplitudes and vice versa. In our 
work, Rcf. [4] , wc demonstrate that all amplitudes of high energy physics are 
some linear combinations of Veneziano amplitudes''^. 

Provided that Eq.(2.4) holds, we shall call such a sequence a partition and 
shall denote it as n =(no, ■■■,nk)- If n is partition of N, then wc shall write n h 
A^. It is well known, e.g. see Ref.s[46,47], that there is one-to-one correspondence 
between the Young diagrams and partitions. We would like to use this fact in 
order to design a partition function associated with Veneziano (and Vcncziano- 
like) amplitudes'^. Clearly, such a partition function should also make physical 
sense. Hence, we would like to provide some qualitative arguments aimed at 
convincing our readers that such a partition function does exist and is physically 
sensible. 

We begin with observation that there is one-to-one correspondence between 

the Young tableaux and directed random walks. Furthermore, it is possible 
to map bijectively such type of random walk back into Young diagram with 
only two rows, e.g. read Ref.[47], page 5. This allows us to make a connection 
with spin chains at once. In this work wc arc not going to use this route 
to spin chains in view of simplicity of the alternative path described in this 
section'^. For this purpose we need to consider a square lattice and to place 
on it the Young diagram associated with some particular partition. Let us 
choose for this purpose some n x m rectangle'"^ so that the Young diagram 
occupies the left part of this rectangle. We choose the upper left vertex of the 
rectangle as the origin of xy coordinate system whose y axis (South direction) 
is directed downwards and x axis is directed Eastwards. Then, the South-East 
boundary of the Young diagram can be interpreted as directed (that is without 
self-intersections) random walk which begins at (0, — to) and ends at (n, 0). 
Evidently, such a walk completely determines the diagram. The walk can be 
described by a sequence of O's and I's. Say, for the x— step move and 1 for the 
y— step move. The totality J\f of Young diagrams which can be placed into such 
a rectangle is in one-to-one correspondence with the number of arrangements 
of O's and I's whose total number is m + n. Recalling the Fermi statistics, the 
number J\f can be easily calculated and is given by TV = {m + n)\/m\n\'^. It 

^For reasons explained in Section 5, it is also called the Poisson-Dirichlet (P-D) distribu- 
tion. 

^"This mathematical result was inspired by the observation that, at least perturbatively, all 
these amplitudes are made of some linear combinations of products of Euler gamma functions 
(with accuracy up to some logarithmic terms). 

'^-'^Tliis task should not be confused with the task of connecting the P-D distributions with 
Veneziano amplitudes to be discussed in Section 5. Alternative pathway through methods of 
asymptotic combinatorics and representation theory [48] will be treated in a separate publi- 
cation. 

^^Nevertheless, this option should not be left underappreciated in view of its immediate 
relevance to Hecke algebra representations, braid groups, etc. e.g. see Ref.s[2] and[49]. 
^•^Parameters n and rh will be specified shortly below. 

^*We have suppressed the tildas for n and m in this expression since these parameters are 
going to be redefined below anyway. 
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can be represented in two equivalent ways: 

(n + l)(n + 2) ■■■{n + m) 



(m + n)!/m!n! 



(m + l)(m + 2) • • • (n + m) 



n! 



n + m 

m 

m + n 
n 



(2.5) 



Let now p(-/V; k, m) be the number of partitions of N into < k nonnegative 
parts, each not larger than m. Consider the generating function of the following 
type 



J^{k,m\q)=Y,PiN;k,m)q'', 



(2.6) 



Af=0 



where the upper limit S will be determined shortly below. It is shown in 

k + m _ k + m 
m ^ k 



Refs. [44,46] that J^{k,m \ q) 
k + m 



where, for instance, 



m 



From this result it should be clear that the expression 

k + m 



k + m 
m 



k + m 
m 



q=l 



is the analog of the binomial coefficient 



In literature [44,46] this 



q— analog is known as the Gaussian coefficient. Explicitly, it is defined as 



a—1 



1) ■ • • (q 



a-b+l 



1) 



(<zf-l)(g^-i-l)-..(rj-l) 



(2.7) 



for some nonegative integers a and b. From this definition we anticipate that the 
sum defining generating function J^{k, m\ q) in Eq.(2.6) should have only finite 
number of terms. Eq.(2.7) allows easy determination of the upper limit S in 
the sum given by Eq.(2.6). It is given by km. This is just the area of the kxm 
rectangle. In view of the definition of p{N; k, m), the number m = N — k. Using 

\ N 

this fact Eq.(2.6) can be rewritten as T{N,k \ q) = 



k 



.This expression 



happens to be the Poincare' polynomial for the Grassmannian Gr{m, k) of the 
complex vector space C^of dimension N as can be seen from page 292 of the 
book by Bott and Tu, Ref.[50]^^. From this (topological) point of view the 
numerical coefficients, i.e. p{N; k,m), in the q expansion of Eq.(2.6) should be 
interpreted as Betti numbers of this Grassmannian. They can be determined 



^^On page 15 of the book by Stanley, Ref.[46], one can find that the number of solutions 

(n + A; — 1 \ 
1^ I j while the number 



of solutions in nonnegative integers to + ... + = n is 



k )■ 



Careful reading of 



Page 15 indicates however that the last number refers to solution in nonnegative integers of 
the equation xo + ■■■ + = n. This fact was used essentially in Eq.(1.21) of Part I. 

^®To make a comparison it is sufficient to repla<;e parameters and n in Bott and Tu book 
by q and N. 
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recursively using the following property of the Gaussian coefRcientsdescribed in 
Ref.[46], page 26, 



n 
k- 



n 
k + 1 



n 
k 



(2.8) 



provided that 



n 




= 1. We refer our readers to Part II for rigorous mathemat- 



ical proof that J^{N, k \ q) is indeed the Poincare' polynomial for the complex 

Grassmannian. With this fact proven, wc notice that, due to relation m = N—k, 
it is sometimes more convenient for us to use the parameters m and k rather 
than N and k. With such a replacement we obtain: 



T{k,m I q) 



k + m 
k 

,m-\-i 



k-\-7n 



k-\-m—l 



1) • • • {q 



m+l 



1) 



M 1-0' 



(2.9) 



This result is of central importance. In our work, Part II, considerably more 
sophisticated mathematical apparatus was used to obtain it (e.g. see Eq.(6.10) 
of this reference and arguments leading to it). 

In the limit q ^ 1 Eq.(2.9) reduces to A/" as required. To make connections 
with results known in physics literature we need to rescale q's in Eq.(2.9), e.g. 
let q = t~ . Substitution of such an expression back into Eq.(2.9) and taking 
the limit t ^ I again reproduces Af in view of Eq.(2.5). This time, however, 
we can accomplish more. By noticing that in Eq.(2.4) the actual value of A'' 
deliberately is not yet fixed and taking into account that m, = N — k, wc can 
fix N by fixing m. Specifically, we would like to choose m = 1 • 2 • 3 • • • fc and 
with such a choice we would like to consider a particular term in the product, 
Eq.(2.9), e.g. 

1 - i^+x 

S{i) = (2.10) 

In view of our "gauge fixing" the ratio m/i is a positive integer by design. This 
means that we are having a geometric progression. Indeed, if we rescale t again 
, e.g. t ^t^, we then obtain 



S{i) = + 



t 



2m 



(2.11) 



with m = Y- Written in such a form the above sum is just the Poincare' poly- 
nomial for the complex projective space CP'". This can be seen by comparing 
pages 177 and 269 of the book by Bott and Tu, Ref.[50]. Hence, at least for some 
m's, the Poincare' polynomial for the Grassmannian in just the product of the 
Poincare' polynomials for the complex projective spaces of known dimensionali- 
ties. For TO just chosen, in the limit t —>■ 1, we reobtain back the number Af as 
required. This physically motivating process of gauge fixing we just described 
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will be replaced by more rigorous mathematical arguments in the rest of this 
paper. Rigorous mathematical arguments causing factorization of the Poincare' 
polynomial can be found, for instance, in Ch-3 of lecture notes by Schwartz, 
Ref.[51]. The relevant physics emerges by noticing that the partition function 
Z[J) for the particle with spin J is given by, e.g. see Ref. [52], 

Z{J) = trie^^"^"^) = e"-^ + e"^'^-^^ + ■■■+ e'"'' 

= 6^=^(1 + e"" + e-^'^ + • • • + e-2=^), (2.12) 

where c is known constant. Evidently, up to a constant, Z{^J) ~ S(i). Since 
mathematically the result, Eq.(2.12), is the Weyl character formula, this fact 
brings the classical group theory into our discussion. More importantly, because 
the partition function for the particle with spin J can be written in the lan- 
guage of N=2 supersymmetric quantum mechanical model^'', as demonstrated 
by Stone, Ref. [52] and others, Ref. [53], the connection between the supersym- 
metry and the classical group theory is evident. It was developed to a some 
extent in Part III. 

In view of arguments presented above, the Poincare' polynomial for the 
Grassmannian can be interpreted as a partition function for some kind of a spin 
chain made of apparently independent spins of various magnitudes'^. These 
qualitative arguments we would like to make more mathematically and physi- 
cally rigorous. The first step towards this goal is made in the next section. 

3 Connection with the Polychronakos-Frahm spin 
chain model 



The Polychronakos-Prahm (P-F) spin chain model was originally proposed by 
Polychronakos and described in detail in Ref. [55]. Frahm [56] motivated by the 
results of Polychronakos made additional progress in elucidating the spectrum 
and thermodynamic properties of this model so that it had become known as 
the P-F model. Subsequently, many other researchers have contributed to our 
understanding of this exactly integrable spin chain model. Since this paper 
is not a review, we shall quote only those works on P-F model which are of 
immediate relevance. 

Following Ref. [55] , we begin with some description of the P-F model. Let 
(7° [a = 1, 2, ...,n^ — 1) be SU{n) spin operator of i-th particle and let the 
operator atj be responsible for a spin exchange between particles i and j, i.e. 

^ij = ^+j:<^r (3.1) 



^^We hope that no confusion is made about the meaning of N in the present case. 

^®In such a context it can be vaguely considered as a variation on the theme of the Polyakov 
rigid string (Grassmann a model, Ref. [54], pages 283-287), except that now it is exactly 
solvable in the qualitative context just described and, below, in mathematically rigorous 
context. 
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In terms of these definitions, the Calogero-type model Hamiltonian can be writ- 
ten as [57,58] 



n = \ UpI + co'x^) + E (3.2) 

^ i i<j (Xi — Xj) 



where I is some parameter. The P-F model is obtained from the above model 
in the limit I — > ±00 . Upon proper rescaling of Ti. in Eq.(3.2), in this limit one 
obtains: 

Hv-^ = -sign{l) E , ,2 , (3.3) 

i<3 [Xi - Xj) 

where the coordinate Xi minimizes the potential for the rescaled Calogero model^^, 
that is 

= E -^—3- (3-4) 

i<d \Xi - Xj) 

It should be noted that H-p^yr is well defined without such a minimization, that 
is for arbitrary real parameters ccj. This fact will be further explained in Section 
4. In the large I limit the spectrum of H is decomposable as 

E = Ec+ lE-p-^, (3.5) 

where Ec is the spectrum of spinless Calogero model while E-p_jr is the spec- 
trum of P-F model. In view of such a decomposition, the partition func- 
tion for the Hamiltonian H at temperature T can be written as a product: 
Z-n{T) =Zc{T)Z'p-j^{T/l). Prom here, one formally obtains the result: 

It implies that the spectrum of P-F spin chain can be obtained if both the total 
and the Calogero partition functions can be calculated. In Ref. [55] Polychron- 
akos argued that ZciT) is essentially a partition function of N noninteracting 
harmonic oscillators. Thus, we obtain: 

Zc{N- T) = n 7^, q = exp(-/?), /? = (fesT)"^ . (3.7) 

Furthermore, the partition function Z-ffiT) according to Polychronakos can be 
obtained using Zc{N\T) as follows. Consider the grand partition function of 
the type 

00 / 00 



H = E Zn{N-T)y'' = E Zc{L;T)y'^ , (3.8) 

Af=0 \L=0 / 

where n is the number of flavors^". Using this definition we obtain: 

n 

Zn{N;T)= \{Zc{ki;T). (3.9) 

Eifei=JVi=l 



i^The Calogero model is obtainable from the Hamiltonian, Eq.(3.2), if one replaces the spin 
exchange operator Cij by 1. Since we are interested in the large I limit, one can replace the 
factor 1(1 — 1) by P in the interaction term. 

^"That is n the same number as n in SU(n). 
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Next, Polychronakos identifies Z„(iV; T) with Z-h(T). Then, with help of Eq.(3.6) 
the partition function Z-p-j^{T) is obtained straightforwardly as 

n (1 - q') 

i=lr=l 

Consider this result for a special case n = 2. It is convenient to evaluate the 
ratio first before calculating the sum. Thus, we obtain: 

N 

n n(i-9'^) 

(3.11) 

where the Poincare' polynomial J-{N, k \ q) for the Grassmanian of the complex 
vector space of dimension N was obtained in the previous section. Indeed, 
Eq.(3.11) can be trivially brought into the same form as that given in our 
Eq.(2.9) using the relation m, + k = N. To bring our Eq.(2.9) in correspondence 
with Eq.(4.1) of Polychronakos, Ref.[55], we use the second line in Eq.(2.9) in 
which we make a substitution: m = N — k. After this replacement, Eq.(3.10) 
acquires the form 



JV k 



4- AN; T) = E n ' ^l.i" (3-12) 



N-i+l 



fe=0 4=0 

coinciding with Eq.(4.1) by Polychronakos. This equation corresponds to the 
ferromagnetic version of the P-F spin chain model. To obtain the antiferromag- 
netic version of the model requires us only to replace qhy in Eq.(3.12) and 
to multiply the whole r.h.s. by some known power of q. Since this factor will 
not affect thermodynamics, following Prahm, Ref.[56], we shall ignore it. As 
result, we obtain 



fe=0 i=0 



in accord with Prahm's Eq.(21). This result is analyzed further in the next 
section. 
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4 Connections with WZNW model and XXX 
s=l/2 Heisenberg antiferromagnetic spin chains 



4.1 General remcirks 



To establish these connections we follow work by Hikami, Ref.[59]. For this 
purpose, we introduce the notation 

n 

(9)n = n(l-5') (4-1) 

i=l 

allowing us to rewrite Eq.(3.13) in the equivalent form 

(4.2) 

Consider now the limiting case {N oo) of the obtained expression. For this 
purpose we need to take into account that 



lim 

AT— >C!0 



N 
k 



To use this asymptotic result in Eq.(4.2) it is convenient to consider separately 

the cases of N being even and odd. For instance, if N is even, we can write: 
= 2m. In such a case we can introduce new summation variables: I = k — m 
and/or I = m — k. Then, in the limit N ^ oo (that is if m— > oo) we obtain 
asymptotically 

Z;/_^(oo;T) = ^ £ /. (4.4a) 
in accord with Ref.[59]. Analogously, if TV = 2m + 1, we obtain instead 

^ oo 

^;^_^(oo;T) = 7-^ E (4.4b) 

According to Mclzcr, Ref. [60], and Kedcm, McCoy and Melzer, Ref. [61], the 
obtained partition functions coincide with the Virasoro characters for SUi (2) 
WZNW model describing the conformal limit of the XXX (s=l/2) antiferro- 
magnetic spin chain, e.g. see Ref. [62]. Even though Eq.s(4.4a) and (4.4b) 
provide the final result, they do not reveal their physical content. This task was 
accomplished in part in the same papers where connection with the excitation 
spectrum of the XXX antiferromagnetic chain was made. Hence, at the physi- 
cal level of rigor the problem of connecting Veneziano amplitudes with physical 
model can be considered as solved. Nevertheless, below we argue that at the 
mathematical level of rigor this is not quite so yet. This conclusion concerns not 
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only problems dicusscd in this paper but, in general, the connection between 
the WZNW models, spin chains and K-Z equations. 

It is true that K-Z equations and WZNW model are inseparable from each 
other as explained, for example, in Ref.[62] but the extent to which spin chains 
can be directly linked to both the WZNW models and K-Z equations still re- 
mains to be discussed. For the sake of space, we shall discuss only the most 
essential facts leaving (with few exceptions) many details and proofs to litera- 
ture. 

Following Varchenko, Ref.[9], we notice that the link between the K-Z equa- 
tions and WZNW models can be made only with help of the Gaudin model, 
while the connection with spin chains can be made only by using the quantum 
version of K-Z equations. Such quantized version of K-Z equations is not im- 
mediately connected with the standard WZNW model as discussed in many 
places, e.g. see Ref.s[9,63]. Therefore, we would like to discuss in some de- 
tail the Gaudin model first and only then its relation to P-F spin chain and, 
accordingly, with the Veneziano model formulated and studied in Part II. 

4.2 Gaudin magnets, K-Z equation and P-F spin chains 

Although theory of the Gaudin magnets plays an important role in topics such 
as Langlands correspondence, Hitchin systems, etc., as explained, for instance, 
in Ref.s[64-66], in this work we do not discuss these topics. Instead, we would 
like to focus only on issues of immediate relevance to this paper. Gaudin came 
up with his magnetic chain model in 1976, Ref.[67], being influenced by earlier 
works of Richardson, Ref.s[68,69] on exact solution of the BCS equations of su- 
perconductivity. This connection with superconductivity will play an important 
role in what follows. 

In physics literature all Gaudin-type models are based on SU{2) algebra of 
spin operators^^. Instead of one Hamiltonian, the set of commuting Hamiltoni- 
ans of the type 

Hi= Y: E<^?^" (4.5) 

is used as discussed in Ref.[71]. In view of the fact that, by construction, 
[Hi,Hj] = 0, and3A^(A^ — 1), the coefficients wfj should satisfy the following 
equations 

<^]k + ^ji^l - = 0- (4-6) 

These equations can be solved by imposing the antisymmetry reqiiircmcnt: 
wfj = —w'ji. It can be satisfied by replacing wfj by the unknown functions 
wfj = /"{zi — Zj) of difference between two new real parameters Zi and Zj. 
It is only natural to make further restrictions based on requirement that the 
2— component of the total spin —'lli '^i is conserved. This causes us to write 

■^^In mathematics literature to be used below [9,63] the SL(2,C) group is used instead of 
its subgroup, SU{2) [70]. 
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wjj = wfj = Xij , and wfj = Yij thus leading to equations 

YijXjk + YkiXjk + XkiXij = 0. (4.7) 
These constraint equations admit the following sets of solutions: 

Xij — Yij = (rational), (4.8a) 

Zi Zj 

^ij ~ ~ — r~ \ ' ^ij ~ cos{zi — Zj) (trigonometric), (4.8b) 

sm [Zi — Zj ) 

Xij = -, Yij = cosh {zi - Zj) (hyperbolic). (4.8c) 

SlUll yZj, '^j) 

While the first solution, Eq.(4.8a), to be used in this work, corresponds to 
the long range analog of the standard XXX spin chain, the remaining two 
solutions correspond to the long range analogs of the XX Z spin chain. 

FoUoving Varchenko, Ref.[9], we are now in the position to write down the K- 
Z equations. For this purpose we combine Eq.s (4.5) and (4.8a) and reintroduce 
the coupling constant g (so that w"j gwfj ) in such a way that the set of K-Z 
equations acquires the form 

d 

{n— - Hi{zi,..., zn))^{zi,..., zn) =0, i = l,...,N, (4.9) 

where K = g~^. This result requires several comments. First, from theory of 
WZNW models it is known that parameter k cannot take arbitrary values. 
For instance, for SUi{2) WZNW model «; = |, e.g. read Ref.[62]. Second, 
we can always rescale z-coordinatcs and to redefine the Hamiltonian to make 
the constant arbitrary small. Apparently, this is asssumed in the asymptotic 
analysis of K-Z equations described in Ref.s[8,9]. Third, if this is the case, 
then such analysis (to be used below) differs essentially from other approaches 
connecting string models with spin chains discussed in the Introduction since 
such a connection was typically made in the limit N oo. Since for SU{N)- 
type models k ~ ^{k + N), in the limit oo we have k ^ oo. The WKB-type 
analysis of K-Z equations of Reshetikhin and Varchenko (to be discussed below) 
fails exactly in this limit. 

With set of K-Z equations defined, we would like now to make a connection 
between the Gaudin and P-F model. To a large extent this was already accom- 
plished in Ref.[72]. Following this reference, we define the spin Calogero (S-C) 
model as follows 

Hs-c = \ J:{p1 + oj'xj) + 9 E 7^r% (4-10) 

i i<j V^* Zj) 

to be compared with H in Eq.(3.2)^^. Using the rational form of the Gaudin 

^^We added the oscillator-type potential absent in the original work, Ref.[72], for the sake of 
additional comparisons, e.g. with Eq.(3.4). In what follows such a constraint is not essential 
and will be ignored. 
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Hamiltonian this result can be equivalently rewritten as 

Hs-c = I J2(p1 + co^xf + i^^\pi,Hi]). (4.11) 
I 

That this is indeed the case can be seen by the following chain of arguments. 

Consider the strong coupling limit {g oo) of Hs-c so that the kinetic 
term is a perturbation. Next, we consider the eigenvalue problem for one of the 
Gaudin's Hamiltonians, i.e. 

if;*W=i;W*W, (4.12) 

and apply the operator ipi to both sides of this equation. Furthermore, consider 

in this limit the combination Hs-c"^^^^- Provided that the eigenvalue problem, 
Eq.(4.12), does have a solution, it is always possible to Fourier expand {ipi'^^'-^) 
using as basis set In such a case we end up with the eigenvalue problem 
for the P-F spin chain in which the eigenfunctions are the same as those for 
the Gaudin's model and the eigenvalues are ipiE^^K Physical significance of this 
result will be discussed in detail below. Before doing so, we have to make a 
connection between the K-Z, Eq.(4.9), and the Gaudin eigenvalue, Eq.(4.12), 
problems. 

Following Ref.s[8,9], we begin by replacing SU{2) spin operators by the 
SL{2,C) = SI2 operators e, / and h obeying commutation relations 

[h,e]=2e; [e,f]=h; [h,f] = -2f. (4.13) 

This Lie algebra was discussed extensively in Part II in connection with design 
of new models reproducing Veneziano amplitudes. In this work, we shall extend 
already obtained results following ideas of Richardson and Varchenko. 

From Ref.[70] it is known that SU{2) is just a subgroup of 5/2- Introduce 
the Casimir element € 5/2 ® s/2 via 

n = e<S:f + f<S:e+^h^h (4.14) 

so that Va; G sh it satisfies the commutation relation [a; 1 + 1 ® x, fJ] = 
inside the U{sl2) ® U{sl2), where (7(5/2) is the universal enveloping algebra of 
5/2- Consider the vector space V = Vi (g) V2 <8) • • • (8> Vjv- An element x € 5/2 
acts on V as follows: a;(8)l(8)l0---(g)l + -- - + l(g)l<g)---(8)x. For indices 
^ ^ i < j < N let f2(''^) : y — > F be an operator which acts as fl on i-th and 
j-th positions and as identity on all others, then the set of K-Z equations can 
be written as 

In the simplest case, this set of equations is defined in the domain U = {{zi, Zn) G 

From now on we shall use Eq.s(4.15) instead of Eq.s(4.9). To connect 
K-Z equations with the XXX Gaudin magnet we shall use the WKB method 
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developed by Reshetikhin and Varchenko, Ref. [8], and summarized in lecture 
notes by Varchenko, Ref. [9]. Following these authors, we shall look for a solution 
of Eq.(4.15) in the form (k ^ 0) 

$(z, k) = e^^(^){/o(z) + «:/i(z) + •••}, (4.16) 

where z = {z\, zn}, S{z) is some scalar function (to be described below) and 
fj{z), j = 0, 1,2, are valued functions. Provided that the function S is 
known, valued functions can be recursively determined (as it is done in the 

WKB analysis). Specifically, given that Hi = ^ , we obtain: 

jj^i Zi — Zj 

Hifoiz) = g^/oiz), (4.17) 
to be compared with Eq.(4.12). Next, we get 

H./,(z) = f/,W + f, (4.18) 

and so on. Since the function S'(z) (the Shapovalov form) plays an important 
role in these calculations, we would like to discuss it in some detail now. 

4.3 The Shapovalov form 

Consider the following auxiliary problem. Let A{x) and B{x) be some pre 
assigned polynomials of degree n and n — 1 respectively. Find a polynomial 
C{x) of degree n — 2 such that the differential equation 

A{x)y" (x) - B{x)y'{x) + C{x)y{x) = (4.19) 

has solution which is polynomial of preassigned degree k. Such polynomial 
solution is called the Lame' function. Stieltjes, Ref.s[8,9], proved the following 

Theorem 4.1. Let A and B he given polynomials of degree n and n — 1, 
respectively so that B{x)/A{x) = Y^^^i ~ - Then there is a polynomial C 

of degree n — 2 and a polynomial solution y(x) = Y['i=i{^ ^ ^i) of Eq.(4.19) if 
and only if x =(xi, ...jXk) is the critical point of the function 

k n 

$fe,„(xi,...,a;fe;2;i,...,2„) = ]J]J(ij- -Zj)"'"" ]J {x^-Xjf. (4.20) 

j=l i=l l<i<j<fe 

Definition 4.2. A point x is critical for $(x) if all its first derivatives vanish 
at it. 
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Wc would like now to make a connection between the Shapovalov form S 
and results just obtained. S is symmetric bilinear form on previously introduced 
space V such that S{v,v) = 1, S{hx,y) = S{x,hy), S{ex,y) = S{x,fy), where 

h, e,f are defined in Eq.(4.13). Furthermore, ^(^(.ti <E)X2),yi 02/2) = S{xi (8> 
X2,^{yi ® 2/2)) Va;i,yi e Vi and Vx2,?;2 £ V2. As result, we obtain: 

S{HiX, y) = S{x, Hiy) Vx, yeV. (4.21) 

Next, let m be some nonnegative integer and Vm be the irreducible Verma 
module with the highest weight m and the highest weight singular vector Vm, 

i. e. 

hvm = mvyn, evm = 0. (4.22) 

Consider a tensor product V = V®^ = Vmi®---®Vmr, so that M = (mi, m„). 
VKn- vectors Vrm^ fvmi,pVmi,- ■ ■,f"^*Vmi form a basis of X^. ^^so that the 
Shapovalov form is orthogonal with respect to such a basis and is decomposable 
as S* = S'mi ® • • • ® -Smn Let, furthermore, J = (ji, j„) be a set of nonnegative 

integers such that ji H \-jn = k, where k is the same as in Eq.(4.20), and < 

ji < This allows us to define the set of vectors /'^fM ~ P^'Vjni®---®f-^'^Vmn 
. These vectors {/"'i'm} are by construction orthogonal with respect to the 
Shapovalov form and provide a basis for the space . Introduce the weight 
of a partition A as \A\ = oi + 02 + ••• then, in view of Eq.(4.22), we define the 
singular vector J^'vm via 

Kf'vM) = {\M\ - 2 \J\)f'vM , e{fvM) = (4.23) 

of weight \M\ — 2 | J| The Bethe ansatz vectors V for the Gaudin model can 
be defined now as 

V(x,z)=^Aj(x,z)/''^;m, (4.24) 

,7 

where x is a critical point of $(x, z) defined by Eq.(4.20) while the function 
Aj{x,t) is defined as follows 

k 

M^,-)= E n^-Z^ (4-25) 



<Ter(k;J) i=l 



with V{k;J) being the set of maps a from the {l,...,fc} to {l,...,n}. Finally, 
using these definitions it is possible to prove that 

5(V(x,z),V(x,z)) = det ( ]n^k,n{xi,-,Xk]Zi,...,Zn)). (4.26) 



■^^ According to Ref. [9] in all subsequent calculations it is sufficient to use the finite Verma 
module, i.e. Lm = Vm/ < f"^^^Vm > ■ This restriction is in accord with our previous 
calculations, e.g. see Part II, Section 8, where such a restriction originates from the Lcfsclictz 
isomorphism theorem used in conjunction with supersymmetric model reproducing Veneziano 
amplitudes. 

^^This fact can be easily undestood from the properties of SI2 Lie algebra representations 
since it is known, e.g see Ref. [9] and Part II, that for the module of highest weight m we have 
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The set of equations determining critical points 

1 d 



$fe,„(xO, zO) dx, 



$fe,„(x(z), z) U=,o=:0 (4.27) 



arc the Bcthc ansatz equations for the Gaudin model. Using these equations the 
eigenvalue, Eq.(4.17), for the Gaudin model now acquires the following form: 

ifi(zO)V(xO,zO) = ^ln$fe,„(x(z), z) |,=,o V(x°,zO). (4.28) 

In the next subsection we shall sudy in some detail the Bethe ansatz Eq.s(4.28).This 
will allow us to obtain eigenvalues in Eq.(4.28) explicitly. 

4.4 Mathematics and physics of Bethe ansatz equations 
for XXX Gaudin model according to Richardson. Con- 
nection with Veneziano model 



Using Eq.(4.20) in (4.28) produces the following set of the Bethe ansatz equa- 
tions: 

n k „ 

E^^=E^— ' i = l.-'fc- (4-29) 

To understand physical meaning of these equations we shall use extensively re- 
sults of two key papers by Richardson, Ref.s [68,69]. To avoid duplications, and 
for the sake of space, our readers are encoraged to read thoroughly these papers. 
Although originally they were written in 60ies having applications to nuclear 
physics in mind, they are no less significant for condensed matter, Ref.[71], and 
atomic physcs, Rcf.[73]. Because of this, only nuclear physics terminology will 
be occasionally used. At the time of writing of these papers, QCD was still in its 
infancy. Accordingly, no attempts were made to apply Richardson's results to 
QCD. Recently, Ovchinnikov, Ref.[74], conjectured that the Richardson-Gaudin 
equations can be useful for development of color superconductivity in QCD. A 
comprehensive review of this topic is given in Ref.[75]. Incidentally, in the same 
paper, Ref. [75], it is emphasized that such type of superconductivity can exist 
only if the number of colors is not too large, e.g. N^. =3. This fact is in accord 
with the remarks made in Section 4.2. regarding the validity of WKB meth- 
ods for K-Z equation in the limit N — > oo. These results clearly favour Abelian 
reduction over ADS/CFT. 

Below, we provide additional mathematically rigorous evidence supporting 't 
Hooft's idea of Abelian reduction of QCD. These results should be considered as 
complementary to that presented in Faddeev's paper, Ref. [27]. For this purpose, 
following Richardson, Ref. [69], we consider the system of interacting bosons 
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described by the (pairing) Hamiltonian^^ 

^ = E,^'"'-fE„'^M;- (4-30) 
Here we have hi = flk , = ''k '^^k ^^'^ M — X] o-kflk- 

k(£k=£i) k(ek=ei) ^(sk=ei) 

It is assumed that the single- particle spectrum {e^} is such that e; < ei+\ VZ 
and that the degeneracy of Z-th level is so that the sums (over k ) each 
contain f2; terms. It is assumed furthermore that the system possesses the 
time-reversal symmetry implying ek = e-k- The operators and ak obey 
usual commutation rules for bosons, i.e. [akja^/] = ^kk'- The sign of the 
coupling constant, in principle, can be both positive and negative. We shall 
work, however, with more physically interesting case of negative coupling (so 
that g in Eq.(4.30) is actually \g\). 

An easy computation using commutation rule for bosons produces the fol- 
lowing results 

[hi,A+]=25u,A+, (4.31a) 
[AuA+]=25ji,{ni + 2h,), (4.31b) 
[hl,A^,] = -25ll,Al. (4.31c) 

If we make a replacement of hi in Eq.s(4.31a) and (4.31.c) by ^ + = ^ 

and keep the same notation in the r.h.s. of Eq. (4.31b) we shall arrive at the 
5^2 Lie algebra isomorphic to that given in Eq.(4.13). The same Lie algebra 
was uncovered and used in our Part II for description of new models describing 
Veneziano amplitudes. Because of this, we would like now to demonstrate that 
the rest of arguments of Part II can be implemented in the present context thus 
making the P-F model (which is derivative of the Richardson-Gaudin XXX 
model) correct model related to Veneziano amplitudes. 

Following Richardson, Ref.[69], we notice that the model described by the 
Hamiltonian, Eq.(4.30), and algebra, Eq.s(4.31), admits two types of excitations: 
those which are associated with the unpaired particles and those with coupled 
pairs. The unpaired z^— particle state is defined by the following two equations 

h\^,>=v\^,>, (4.32) 

Ai I VL (4.33) 

Here, h = ni so that, in fact, 

fii I >= vi\^^> (4.34) 

and, therefore, v = ui . Furthermore, 

^^In the paper with Sherman, Ref.[68], Richardson explains in detail how one can map the 
fermionic (pairing) system into bosonic. 



20 



Following Richardson, wc want to demonstrate that parameters si in Eq.(4.35) 
can be identified with parameters zi in Bethe Eq.s (4.29). Because of this, the 
eigenvalues for the P-F chain are obtained as described in Section 4.2., that is 

^r^' = |-E,-'-'=-^ (4-36) 

These are the eigenvalues of fii defined in Eq.(4.34). Furthermore, this eigen- 
value equation is exactly the same as was used in Part II, Section 8, with purpose 
of reproducing Veneziano amplitudes. Moreover, Eq.s(4.32) and (4.33) have the 
same mathematical meaning as Eq.s (4.23) defining the Vcrnia module. Because 
of this, we follow Richardson's paper to describe this module in physical terms. 
By doing so additional comparisons will be made between the results of Part II 
and works of Richardson. Since the Hamiltonian, Eq.(4.30), describes two kinds 
of particles: a) pairs of particles (whose total linear and angular momentum is 
zero) and, b) unpaired particles (that is single particles which do not interact 
with just described pairs), the total number of (quasi) particles is n = A'' + u'^^. 

Since we redefined the number operator as ^ +n; = — = N(, we expect that , 

^ 4 

once the correct state vector describing excitations is found, Eq.(4.30) should be 
replaced by the analogous equation for N( whose eigenvalues will be ^ + ^(.'^^ 
A simple minded way of creating such a state is by constructing the following 
state vector j4j^ •••A^ \(p^>. This vector does not possess the needed symmetry 
of the problem. To create the state vector (actually, the Bethe vector of the 
type given by Eq.(4.24)) of correct symmetry one should introduce a linear 
combination of operators according to the following prescription: 

B+=J2ua{l)A+, a = l,...,N (4.37) 

with constants Ua(l) to be determined below. The (unnormalizcd) Bethc-type 
vectors are given then as \ -ip >= ■ ■ ■ \ > and, accordingly, instead of 
Eq.(4.35), we obtain 

^ I ^ >= ^'^') I ^ > +[^' Bt ■■■B+]\p,> . (4.38) 

The task now lies in calculating the commutator and to determine the constants 
Ua{l). Details can be found in Richardson's paper, Ref.[69]. The final result looks 
as follows 

H\i}> -E\il)> (4.39) 

N 

= ^(n^7)I3^^[(2£;-^a)Wa(0 + 13(^i'+2n,0Wa(O+4<? ^ M^„] | (/P, > . 
a=l 75^a I I' PiPjtot) 

^®In Richardson's paper we find instead: n = 2N + v. This is, most likely, a misprint as 
explained in the text. 

^^These amendments are not present in Richardson's paper but they are in accord with its 
content. 
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By requiring the r.h.s. of this equation to be zero, we arrive at the eigenvalue 
equation 

JV 

I V >= -B I V' >, where E = + ^ (4.40) 

I a=l 

Furthermore, this requirement after several manipulations leads us to the Bethe 
ansatz equations^* 

as well to the explicit form of coefficients Uq,(Z) : Ua{l) = l/(2e; — Ea) and that 
for the matrix elements Mq_^ (since, by construction, Ua{l)u(3{l) = Majyu,a(l) + 
Mfj_aU(j{l)). In the limit g — > we expect E^ —>■ '2si and Jl; — > in accord with 
Eq.s(4.32)-(4.34). Therefore, we conclude that ^ + 1^1 is an eigenvalue of the 
operator N; acting on | ^ > in accord with remarks made before. In the opposite 
limit: g( — > 00, the system of Eq.s(4.41a) will coincide with Eq.(4.29) upon 
obvious identifications: ^ E^, 2ei ^ zi, N ^ k, L ^ n and Qi/2 + 1^1^ mi. 
Next, in view of Eq.s(4.36) and (4.40) we obtain the following result for the 
occupation numbers: 

A JV 

(4.42) 



1 dsi 



Based on the results just obtained, it should be clear that, actually, ' 

JV 

=i'i + ^ so that ^ = J2 Richardson, Ref.[69], cleverly demonstrated 



that the combination ^ must be an integer. 

a=l ^' 

Consider now a special case: A'^ = 1. Evidently, for this case, the derivative 
should also be an integer. For different e-s these may, in general, be 
different integers. This fact has some physical significance to be explained below. 

To simplify matters, by analogy with theory of superconducting grains, 
Ref.[71], we assume that the energy Si can be written as £, = d{2i — L — 1), 
i = 1,2, ...,!/. The adjustable parameter d measures the level spacing for the 
unpaired particles in the limit g ^ 0. With such simplification, we obtain the 
following BCS-type equation using Eq.s(4.41a) (for N = 1): 



i-i2e,-E G' 



(4.43) 



^^It should be noted that in the original paper, Ref. [69], the sign in front of the 3rd term in 
the l.h.s. is positive. This is because Richardson treats both positive and negative couplings 
simultaneously. Equation (4.41a) is in agreement with (3.24) of Richardson-Sherman paper, 
Ref. [68] , where the case of negative coupling (pairing) is treated. 
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Figure 1: Graphical solution of the Eq.(4.44) 



where G is the rescaled coupling constant. Such an equation was discussed in 
the seminal paper by Cooper, Ref. [76], which paved a way to the BCS theory 

L 

of superconductivity. To solve this equation, let now F{E) = ^ ^i{2ei — E)^^ 

1=1 

so that Eq.(4.43) is reduced to 



F{E) = G- 



(4.44) 



This equation can be solved graphically as depicted below, in Fig.l. 

As can be seen from Fig.l, solutions to this equation for G = oo can be 
read off from the x axis. In addition, if needed, for any > 1 the system 
of Eq.s (4.41a) can be rewritten in a similar BCS-like form if we introduce the 
renormalized coupling constant Ga via 



^ 1 



(4.45) 



so that now we obtain: 



F{E^ 



G„ 



(4.41b) 



This system of equations can be solved iteratively, beginning with Eq.(4.44). 
There is, however, better way of obtaing these solutions. In view of Eq.s(4.19), 
(4.20) and (4.27) solutions {Ea} of Eq.(4.41.b) are the roots of the Lamc'-type 
function which is obtained as solution of Eq.(4.19). Surprisingly, this fact known 
to mathematicians for a long time has been recognized in nuclear physics liter- 
ature only very recently, e.g. read Ref. [77]. 
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4.5 Emergence of the Veneziano-like amplitudes as con- 
sistency condition for = 1 solutions of K-Z equa- 
tions. Recovery of pion-pion scattering amplitude 

Since results for the Richardson-Gaudin (R-G) model are obtainable from the 
corresponding solutions of K-Z equations, in this subsection we would like to 
explain why iV = 1 solution of the Bethe-Richardon equations can be linked with 
the Veneziano-like amplitudes describing the pion-pion scattering. In doing so, 
we shall by pass the P-F model since, anyway, it is obtainable from the R-G 
model. 

Thus, wc begin again with Eq.s (4. 14), (4. 15). We would like to look at 
special class of solutions of Eq.(4.15) for which the parameter \ J\ in the Verma 
module, Eq.(4.23), is equal to one. This corresponds exactly to the case = 1. 
FoUoving Varchenko, Ref.[9], by analogy with Eq.(4.20) we introduce function 
$(z, f) via 

$(z,t)= H {z,~Zj) K Y[{t-zif~. (4.46) 

l<i<i<L 1=1 

It is a multivalued function at the points of its singularities, i.e. at the points 
Zi,...,Zl. Using this function, we define the set of 1-forms via 

a;j=$(z,i)-^, j = l,...,L, (4.47) 

t Zj 

and the vector I^'''^ of integrals I^'''^ = (Ii, = (Xy '^i; •••) I^^l) with 7 being 
a particular Pochhammer countour: a double loop winding around any two 
points Za, Zfj taken from the set zi, z^. Deatails can be foimd in Rcf.s[9,63]. 

Wc want now to design the singular Verma module for the K-Z equations 
using Eq.(4.23) and results just presented. Taking into account the following 
known relations: 

a) ef'^Vm = k{m -k + l)/''"^^™, and b) hf'^Vm = (m - 2k)f''vm 

for the Lie algebra SI2, also used in Part II, Section 8, and taking into account 
that in the present {N = 1) case the basis vectors /"^Vm = f^^ ®---® f-'"Vm„ 
acquires the form: f^VM = v„i^ ® • • • ^ fvm,, (E) ■ ■ ■ (E) Vm^, -s = 1, L, provided 
that m-s are the same as in Eq.(4.29) (or (4.46)), the singular vector for such 
a Verma module is given by 

L 

^(7) = '^l8Vmi^ - ■ - ^ fVms (8) • • • (8) Vm„ • (4.48) 

s=l 

In view of the Lie algebra relations just introduced, we obtain e • w = or, 
explicitly, 

L 

J2msls=0. (4.49) 

s=l 
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Hence, for a fixed Pochhammer contour 7 there are L — f independent basis 
vectors {w^}. They represent L — 1 independent solutions of K-Z equation of 
the type k = 1 (oi N = 1). Let now z'^s be ordered in such a way that zi < 
■■■ < zl- Furthermore, in view their physical interpretation described in previous 
section, these z[s can be chosen to be equidistant. Consider then a special set 
of Pochhamer contours {7,} around points Zi and -Jj+i, i = 1,2,...,L — 1, 
and consider the matrix M made of integrals of the type M] = f uij 

then, any (fc = 1)— type solution ^^{i = 1, 2, i — 1) of K-Z equation can be 
represented as 



From linear algebra it is known that in order for these K-Z solutions to be 
independent we have to require that det M ^ 0. The proof of this fact is given 
in the Appendix. Calculation of the determinant of M is described in detail in 
Ref.[9] so that we quote the result: 



where ±A is some known constant^^ and r(a;) is Euler's gamma function. For 
L = 2 without loss of generality one can choose zi = and Z2 = 1, then in the 
determinant thus obtained one easily can recognize the Veneziano-type tt+tt" 
scattering amplitude used in the work by Lovelace, Ref.[78]. We discussed this 
amplitude previously in connection with mirror symmetry issues in our work, 
Ref.[79]. This time, however, we would like to discuss other topics. 

In particular, we notice first that all mesons are made of two quarks. Specif- 
ically, we have: ud for ir^^du for 7r~ and dd for tt*^. These are very much like 
the Cooper pairs with qq quark pairs contributing to the Bose condensate cre- 
ated as result of spontaneous chiral symmetry breaking. As in the case of more 
familiar Bose condensate, in addition to the ground state we expect to have 
a tower of excited states made of such quark pairs. Experimentally, these are 
interpreted as more massive mesons. Such excitations are ordered by their ener- 
gies, angular momentum and, perhaps, by other quantum numbers which can be 
taken into account if needed. Color confinement postulate makes such a tower 
infinite. Evidently, the Richardson-Gaudin model fits ideally this qualitative 
picture. Eq.(4.44) describes excitations of such Cooper- like pairs (even in the 
limit: G — > 00) as can be seen from Fig. 1. In the P-F model the factor plays 
effectively the role of energy as already discussed in this work and in Part II. 
Therefore, in view of Eq.(4.42), it is appropriate to write: fli = f{Ei), with Ei 
being the R-G energies. Although the explicit form of such /—dependence may 
be difficult to obtain, for our purposes it is sufficient only to know that such a 




(4.50) 



detM = ±A 



r(i - SL) . . . r(i - 2^) 
r(i-^) 



(4.51) 



±A= n (2<-2j) « 



l<i,j<L 
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dependence does exist. This then allows us to make an identification: Oj — - 

K 

consistent with Varchenko's results, e.g. compare his Theorem 3.3.5 (page 35) 
with Theorem 6.3.2. (page 90) of Ref.[9]. But, we already established that ^li is 
an integer, therefore, — - should be also an integer. This creates some apparent 

K 

problems. For instance, when \M\ = k, the determinant, detM, becomes zero 
implying that solutions of K-Z equation become interdependent. This fact has 
physical significance to be discussed below and in Section 5. To do so we use 
some results from our Part I. In particular, a comparison between 



-(!))(! H-(^)) (4.52) 



Sm TTZ = TTZ 

fc=l 

and 

^"''^IKi+It))^"^' (4-53) 



ze 



z 



Viz) 
^ ' fc=i 

where C is some known constant, tells us immediately that not only |Af| = 
K will cause detM=0 but also \M\ = K{k + l),k — 0,1,2,... Accordingly, 

the numerator in Eq.(4.51) will create poles whenever — = 1. Existence of 

K 

independent K-Z solutions is not destroyed if, indeed, such poles do occur. 

TTl ' 

These facts allow us to relabel — - as a{s) (or a{t) or a{u), etc.) as it is done in 

K 

high energy physics with continuous parameters s,t,u, ...replacing discrete z's, 
different for different F functions in the numerator of Eq.(4.51). In the simplest 
case, this allows us to reduce the determinant in Eq.(4.51) to the form used by 
Lovelace, i.e. 

detM = -Alil^«i^. (4.54) 
r(l - a{s) - a{t)) ^ ' 

If, as usual, we parametrize a(s) = a(0) + a's, then equation 1 = a(s) + a{t) 
causes the detM to vanish. This also fixes the parameter Q!(0): Q!(0) = 1/2. 
This result was obtained by Adler long before sting theory emerged and is 
known as Adler's selfconsistency condition, Ref.[80]. With such "gauge fixing", 
one can fix the slope a' as well if one notices that the experimental data allow 
us to make a choice: 1 = a(m^).This leads to: o! = ~ 0.885{Gev~^), in 
accord with observations. 

The obtained results are not Hmited to study of excitations of just one "su- 
peconducting" pair of quarks. In princile, any finite amount of such pairs can 
be studied, e.g. see Ref.[81]. In such a case the result for detM is expected to 
become considerably more complicated but connections with one dimensional 
magnets still remain unchanged. We plan to discuss these issues in future pub- 
lications. 
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5 Random fragmentation and coagulation pro- 
cesses the Poisson-Dirichlet distribution and 
Veneziano amplitudes 

5.1 General facts about the Poisson-Dirichlet distribution 

In the Introduction, following Heisenberg, we posed a question: Is combina- 
torics of observational data sufficient for recovery of underlying unique micro- 
scopic model? That is, can we have complete understanding of such a model 
based on information provided by combinatorics? As we demonstrated, espe- 
cially in Section 4 and in Rcf. [37], this task is impossible to accomplish with- 
out imposing additional constraints which, normally, are not dictated by the 
combinatorics only. Even accounting for such constraints, the obtained results 
could be in conflict with rigorous mathematics and physical reality. Last but 
not the least, since Veneziano amplitudes gave birth to string theory one can 
pose another question: Is these Veneziano (or Veneziano- like) amplitudes, per- 
haps corrected to account for particles with spin, contain enough information 
(analytical, number-theoretic, combinatorial, etc.) allowing restoration of the 
underlying microscopic model uniquely? In the most general case the answer 
is: No! This happens in spite of the fact that all amplitudes of high energy 
physics can be made out of linear combination of Veneziano amplitudes (up to 
logarithmic corrections) as discussed in our recent work, Ref.[4]. In the rest of 
this section we explain why this is so. 

We begin with recalling some known auxiliary facts from the probability 
theory. For instance, we recall that the stationary Maxwell distribution for ve- 
locities of particles in the gas is of Gaussian-type. It can be obtained as a 
stationary solution of the Boltzmann's dynamical equation maximizing Boltz- 
mann's entropy^°. The question arises: Is it possible to find (discrete or contin- 
uous) dynamical equations which will provide known probability distributions 
as stable stationary solutions? This task will involve finding of dynamical equa- 
tions along with the corresponding Boltzmann-like entropies which will reach 
their maxima at respective equilibria for these dynamical equations. We are 
certainly not in the position in this work to discuss this problem in full general- 
ity. Instead, we focus our attention only on processes described by the so called 
Dirichlct distributions. These originate from the integral ( e.g. see Eq.(2.8) of 
Part I) attributed to Dirichlct. It is given by 

V{xi,...,Xn+l)= J ■ J ur"^---<""^(l-Ul Unf"+-'-'^dUi---du„. 

Ui>0,..., Un>0 
Ui-\ \-Un<l 

(5.1) 

n 

A random vector (Xi,...,X„) G R" such that Xj > and Yl =1 is 

i=l 

said to be Dirichlct distributed with parameters {xi, ...,a;„;a;„+i), e.g. see Ref. 
^"As discussed in our work, Ref.[82], on the Poinccire' and geometrization conjectures. 
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[83], if the probability density function for (Xi, X„) is given by 

-("'•■■■•"»> = ?g.).:.rU:) <-' ■ ■ • - ■ 

r(a;i H \-Xn+i) J, _i x„.i-i -J J XI, X 

= r(xi) ■ ■ ■ r(a: +i) , provided that M„+i 

= 1 - Ul Un- (5.2) 

From these results it follows that Veneziano condition, Eq.(2.4), and Veneziano 
amplitudes are inseparable from each other. Since Veneziano condition is just 
restatement of energy-momentum conservation, such a requirement should be 
applicable to whatever amplitude of high energy physics. Not surprisingly, there- 
fore, in Ref . [4] , it is demonstrated that this is indeed the case. 

It is of interest to mention other uses of Dirichlct distributions beyond that 
in high energy physics. For this purpose, to get a feeling of just defined dis- 
tribution, we notice the following peculiar aspects of this distribution. For any 
discrete distribution, wc know that the probability pi must be normalized, that 
is X^jPi = 1. Thus, the Dirichlet distribution is dealing with averaging of the 
probabilities! Or, better, is dealing with the problem of effectively selecting 
the most optimal probability. The most primitive of these probabilities is the 
binomial probability given by 

Pm= ( ^ )p'"(l-pr-", m = 0,l,2,....,n. (5.3) 

If X is random variable obeying this law of probability then, the expectation 
E{X) is calculated as 

n 

E{X) = ^ mpm = np = iJ.. (5.4) 

m— 1 

Consider such a distribution in the limit: n ^ oo. In this limit, if we write 
p = ji/n , then the Poisson distribution is obtained as 

, TO 

Pm = ^e-^ (5.5) 
ml 

Next, we notice that m! = r(m+l), furthermore, we replace m by the real valued 
variable a and fihy x. This allows us to introduce the gamma distribution with 
exponent a whose probability density is 

Px{x) = TT^^-'e-- (5.6) 

for some gamma distributed random variable X. Finally, we would like to 
demonstrate how the Dirichlet distribution can be represented through gamma 
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distributions. Since the gamma distribution originates from the Poisson distri- 
bution, sometimes in literature the Dirichlet distribution is cahed the Poisson- 
Dirichlet (P-D) distribution, Ref.[45]. To demonstrate connection between the 
Dirichlet and gamma distributions is relatively easy. Following Kingman, Ref.[45], 
consider a set of positive independent gamma distributed random variables: 
Yi, Yn^i with exponents ai, a„+i. Furthermore, consider Y = Yi + ■ ■ ■ + 
Yn+i and construct a vector u with components: Ui = Then, since '^'^^i Ui 
= 1, the components of this vector are Dirichlet distributed and, in fact, inde- 
pendent of Y. Details of the proof are based on results already discussed in 
Part I and are given in Appendix. 

Such described Dirichlet distribution is an equilibrium measure in various 
fields ranging from spin glasses to computer science, from linguistics to ge- 
netics, from forensic science to economics, etc. Many useful references involv- 
ing these and other applications can be found in Ref.s[38-40]. Furtheremore, 
most of fragmentation and coagulation processes involve the P-D distribution 
as their equilibrium measure. Some applications of general theory of these pro- 
cesses to to nuclear and particle physics were initiated in a series of papers by 
Mekjian, e.g. see Ref.[41] and references therein. Alternative approach to the 
fragmentation-coagulation processes in high energy physics was developed by 
Andersson, Ref.[42], and is known as the Lund model. As results of our re- 
cent work, Ref.[4], indicate, the results of Mekjian, Ref.[41], and that presented 
in this work are already fully compatible with general theory of coagulation- 
fragmentation processes discussed in Ref.s [38-40]. The interconnections between 
the Lund model and general theory of coagulation-fragmentation processes re- 
mains to be investigated. In the meantime, we woiild like to connect general 
results presented in this subsection with those of Section 4. This is accomplished 
below. 

5.2 Quantum mechanics, hypergeometric functions and 
P-D distribution 

In Ref.[2,3] we provided detailed explanation of the fact that all exactly solv- 
able 2-body quantum mechanical problems involve different kinds of special 
functions obtainable from the Gauss hypergeometric funcftion whose integral 
representation is given by 

1 

^^'"'^'"'^^ = r(b)r(c- b) J ^'"'^^ " - ^^)""'^*- (5.7) 



As is well known, the disctete spectrum of all exactly solvable quantum me- 
chanical problems can be obtained only if one can find an appropriate set of 
orthogonal polynomials related to this spectrum. Since all these orthogonal 
polynomials are obtainable from Gauss hypergeometric function, the question 
arises: Under what conditions on coefficients (a, b and c) can infinite hypergeo- 
metric series (whose integral representation is given by Eq.(5.7)) be reduced 
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to a finite (orthogonal) polynomial? This happens, for instance, if we im- 
pose the quantization condition: —a = 0,1,2,... In such a case we can write 
(1 - zt)-"- = Ej^"i(i'")(-l)'(-2*)' and use this finite expansion in Eq.(5.7). In 
view of Eq.(5.2), wc obtain the convergent generating function for the Dirich- 
let distribution. Hence, all known quantum mechanical problems involving 
discrete spectrum are effectively examples of the P-D stochasic processes. Fur- 
thermore, from this point of view quantum mechanics ( also, quantum field 
theory and string theory, e.g. see Ref.[4]) becomes just an applied theory of the 
P-D stochastic processes. For hypergeometric functions of multiple arguments 
this was demonstrated in Ref.[84] only quite recently. Other arguments in favour 
of such an interpretation are developed in Ref.[4]. 

Next, we are still interested in the following. Given these observations, can 
we include the detcrminantal formula, Eq.(4.51), into emerging quantization 
scheme? Very fortunately, this can be done as explained in the next subsection.. 

5.3 Hypergeometric functions, Kummer series and Veneziano 
amplitudes 

In view of just introduced new quantization condition, the question arises: Is 
this the only condition reducing the hypergeometric function to a polynomial 
? More broadly: what conditions on coefficients a, b and c should be imposed 
so that the function F{a, b, c; z) becomes a polynomial? The answer to this 
question was provided by Kummer in the first half of 19th century as discussed 
in Ref. [85] . Incidentally, in the case of K-Z equations such a problem was solved 
only in 2007 in Ref. [86] . We would like to summarize Rummer's results and to 
connect them with the dcterniinantal formula, Eq.(4.51). 

According to general theory of hypergeometric equations of one variable 
discussed in Ref. [85] , the infinite series for hypergeometric function degenerates 
to a polynomial if one of the numbers 

a,b,c—aoT c — b (5.8) 

is an integer. This condition is equivalent to the condition that, at least one of 
eight numbers ±(c — 1) ± (a — 6) ± (a + 6 — c). is an odd number. According 
to general theory of hypergeometric functions of multiple arguments discussed 
in Section 4, the k = 1-type solutions can be obtained using 1-forms given by 
Eq.(4.47) accoimting for a singular module constraint, Eq.(4.49), in the form 
given by Eq.(4.46). In the case of Gauss- type hypergeometric functions, rela- 
tions of the type given by Eq.(4.49) were known already to Kummer. He found 
24 interdependent solutions. Evidently, this number is determined by the num- 
ber of independent Pochhamer contours as explained in Ref.s[9, 85]. Therefore, 
among these he singled out 6 (generating these 24) and among these 6 he estab- 
lished that every 3 of them are related to each other via equation of the type 
given by Eq.(4.49). 

We denote these 6 functions respectively as ui, ue- Then, we can represent, 
say, U2 and uq using ui and U5 as the basis set. We can do the same with ui 
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and U5 by representing them through U2 and uq and, finally, we can connect M3 
and U4 with ui and M5. Hence, it is sufficient to consider, say, U2 and uq. Thus, 
we obtain: 

W2 A / Ml Ml \fui 
ue J [ Ml M| j I M5 



(5.9) 



with Ml = na + b-c+mi~c) ^ r(a + 6+l-c)r(c-l) ^ 

1 r(a + i-c)r(6-c + i)' 2 r(a)r(6) ' ^ 

r(c + i-a-fe)r(i-c) r(c + i-a-5)r(c-i) ^, , . 

— — ^- , / \, '-■,Mi = ^ , / '-. The determinant of 

r(l-a)r(l-6) ' ^ T{c-a)T{c-b) 

this matrix becomes zero if either two rows or two columns become the same. 
For instance, we obtain: 

r(a)r(b) _ r(a - c + i)r(b - c + 1) r(c - Q)r(c - 5) _ r(i - a)r(i - h) 
r(c-i)~ r(i-c) r(c-i) ~ r(i-c) 

(5.10) 

The condition c = 1 in Eq.(5.10) causes two solutions of hypergcomctric 
equation to degenerate into one polynomial solution in accord with general 
theory. 

Appendix 

A. Linear independence of solutions of K-Z equation 

Linear independence of solutions of K-Z equation is based on the following 
arguments. Consider change of the basis 

e-'=A^e\*,.j = l,2,...,n (A.l) 

in R". Using this result, consider the exterior product 

eiA---Ae" = [detA]e^ A---Ae". (A.2) 

Next, suppose, that the vectors are lineraly-dependent. In particular, this 
means that 

e" = aigi + • • • + a„_ie"-i (A.3) 
for some nonzero a[s. Using this expansion in Eq.(A.2) we obtain 

e^ A • • • A e"-i A (aie^ + • • • + a„_ie"-i) = 0, (A.4) 

implying [det A] =0. Convesely, if [det A] 7^0 then, vectors are linearly in- 
dependent. 



B. Connections between the gamma and Dirichlet distributions 

Using results of our Part I (especially Eq.(3.27)), such a connection can be 
easily established. Indeed, consider n + 1 independently distributed random 
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gamma variables with exponents ai, Q:„+i. The joint probabiHty density for 
such variables is given by 

Pr.,...r.,M, .„.r) = ^ • • • J^^^"^ • • • s^V-'- (B-D 

Let now Sj = tit, where ti are chosen in such a way that YZtlU = 1- Then, 
using such a substitution in Eq.(B.l), we obtain at once: 

oo 

p^^,...^^^^{tu...,tn+i) = [ft"-'e-']-^----^ (B.2) 



Since a = ai + • • • + a„+i, we also obtain: J t"~^e~* — ^{oti + • • • + Un+i), 



implying that the density of probability given by Eq.(B.2) is indeed of Dirichlet- 
type given by Eq.(5.2) of the main text. 
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